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Abstract 

We calculate the two-loop electroweak corrections to the fermionic 
decay width of the Higgs boson at leading order in ran- The sum 
of one-loop and two-loop radiative corrections turns out to be at 6% 
level over the whole range of mn where the perturbation theory is sup- 
posed to be a sensible approximation. We show that the perturbative 
approach breaks down at ~ 1.4 TeV, and address the question of 
its relevance beyond this limit. 



That the selfinteraction of the electroweak symmetry breaking sector of 
the standard model increases with the mass of the Higgs boson, and that 
this eventually leads to the breakdown of the perturbation theory, belongs 
to the common wisdom. To establish the limits of perturbation theory is, 
however, a more delicate matter. 

The coupling constant of the Higgs sector, £^^ L , becomes of order unity 
for mu ~ 1-5 TeV. This is roughly where one expects the problems related 
to the asymptotic nature of the perturbative series to show up. 

The perturbative solution for the S matrix is formally unitary, but its 
approximations of any finite order are not. The unitarity violation effects 
are of higher order in the coupling constant. If the perturbation theory 
fails to converge satisfactorily, they become numerically large, impairing 
the relevance of the perturbative amplitudes. The well-known partial wave 
analyses of longitudinal vector boson scattering at tree level by Dicus and 
Mathur Q , and B.W. Lee, C. Quigg and H. Thacker || show that the the 
unitarity bounds are exceeded for mu ~1 TeV. Durand, Johnson and Lopez 

U extended this analysis to include the one-loop radiative corrections, 
and extracted a bound of ~ 400 GeV. This much smaller value is, however, 
to be traced back to their criterion for unitarity violation, which is stronger 
than that of refs. |2) , rather than to genuine one-loop effects. 

It is, of course, hard to draw out unambiguous conclusions about the 
breakdown of the perturbative series if only its first two terms are known. 
The first nontrivial criterion is to compare the magnitude of the one-loop 
and two-loop radiative corrections, and to establish where the series starts 
to diverge. 

Here one must distinguish between radiative corrections to low energy 
parameters, and to processes involving the symmetry breaking scalars. The 
latter are expected to be larger, since the contributions which are of leading 
order in the Higgs mass cancel in the radiative corrections to low energy 
parameters 0, |j. For instance, the two-loop radiative corrections to the p 
parameter and to the vector boson selfcouplings are quadratic in m#, up to 
logarithms of mu- They become comparable to the one-loop logarithmic 
corrections only for Higgs bosons as heavy as 3 — 4 TeV || 0]. 

On the other hand, performing multiloop radiative corrections to pro- 
cesses which involve the Higgs boson is a difficult task. Already in the case 
of the two-loop massive diagrams no general analytical solution exists for 
nonvanishing external momenta. One is therefore forced to calculate such 
diagrams at least partly numerically. 

We calculate in this paper the leading two-loop corrections to the fermio- 
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nic width of the Higgs boson, which grow quartically with mn- We use the 
method described in ref. |J to calculate the massive two-loop diagrams. 
The two-loop correction becomes larger than the one-loop one if mn > 1.4 
TeV, signaling the breakdown of the perturbation theory. This is compatible 
with the results of ref. ||, where the two-loop corrections to the shape of 
the Higgs resonance were calculated and shown to become large for mu ~ 
1.2 TeV. This is also roughly the value one would expect by considering 
the magnitude of the coupling constant of the Higgs sector. On the other 
hand, this disagrees with the conclusions of ref. about the breakdown 
of perturbation theory at 380 GeV []. In view of this disagreement, special 
emphasis is put on checks of both the analytical and numerical parts of the 
calculation. 

Since we are interested in the leading m# effects, the natural choice 
is to work in Landau gauge. Power counting arguments show that in this 
gauge only the Feynman diagrams with scalar propagators survive || o[ . 
The diagrams with fermions, gauge bosons, and Fadeev-Popov fields need 
not be taken into account since they are not of leading order. However, we 
effectively evaluate the diagrams in a nearly-Landau gauge, that is, we keep 
a small, finite gauge parameter £ throughout the calculation, and take the 
limit £ — > in the final result. This is necessary at two-loop level in order 
to avoid the arbitrariness of J d n p ^ within the dimensional regularization. 
One can convince oneself that this is a regular limit, in the sense that no 
finite contributions from the diagrams with gauge or Fadeev-Popov ghosts 
survive. This procedure provides also a good check of the calculation, since 
all £ dependence must cancel in the final results. 

Within this framework, the leading corrections to the Yukawa couplings 
— hig Hff are a pure renormalization effect. They are simply given by 

a factor ( — ^ ) 2 , so one only needs to evaluate selfenergy diagrams. 
1 — 

Denoting the bare Higgs and Goldstone fields by Ho, and z°, one 
writes down the Higgs sector of the standard model as: 

1 The two-loop radiative corrections of ref. |u| are one order of magnitude larger, which 
is also roughly the order of magnitude of the cancellations between different diagrams. Ref. 
Jic| ] probably misses these cancellations because of a few diagrams evaluated incorrectly, 
two of which were already identified in ref. ||. One also remarks that their dSo, which 
they identify with zero, is a quantity which diverges logarithmically as the mass of the 
Goldstone bosons tends to zero. The proof that dSo = given in ref. is based 

on illegally splitting this integral into two contributions whose definition domains do not 
overlap. 



2 



C = ^(^F )(^iIo) + ^(Vo)(^o) + (^0(^ -) 

2 m Ho 1 , + - , 1 2 , 1 ^2 | 2m W „ , 4Jt 2 
~g Z^o [< W + ^0 + + — — #0 + m2 ] ,(1) 

m W H 2 2 g 2™hq_ 

Wq 

where 5t is the tadpole counterterm needed to ensure that the vacuum ex- 
pectation value of the Higgs field does not receive quantum corrections. It 
is related to the scalar selfenergy at zero momentum. One further splits the 
Lagrangian of eq. 1 into a renormalized Lagrangian and counterterms: 

tj nr^/2 yj 

HQ — Zj H n 

rrl/2 

Zo = Zq z 

ryl/2 

Wq = Z(j W 

2 2 x 2 

m Ho = m H -Sm H 

m w = m w ~ ^ m w i ( 2 ) 

and fixes the counterterms by using physical, on-shell renormalization con- 
ditions with field renormalization. Note that the gauge coupling constant 
g does not get renormalized at leading order in the coupling constant of 

2 

the Higgs system g 2 —^. One can express g in terms of the Fermi coupling 

m w 

constant and the mass of the charged vector boson as g 2 = 4^/2ra^ Gf, 
with G F = 1.16637 • 1(T 5 GeV~ 2 , and m w = 80.6 GeV. 

The necessary one-loop counterterms at 0(e) are given in ref. ||]. 

One can now proceed with the actual two-loop calculation. The main 
task is to evaluate the W boson selfenergy at zero momentum for extracting 
the W mass counterterm, and the on-shell derivative of the Higgs selfenergy 
for calculating the field renormalization constant of the Higgs boson. This 
involves the selfenergy topologies of. fig. 1. An algebraic computer program 
generates all relevant Feynman diagrams, and uses the techniques of ref. ||] 
to calculate them. In the case of the selfenergy of the W boson it is possi- 
ble to obtain an analytic result because the external momentum vanishes. 
For the wave function renormalization of the Higgs boson, some numerical 
integrations have to be performed. The results are: 
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r 2(2-loop) _ . 2rn 2 H 2 m 2 H m 2 w 3 1 J_ 3 
19 m 2 v j W J (16vr2)2 [ 32e 128 + 32 7 

vr 2 3\/37r _ 3\/3 7T 
192 64 16 4 jJ 



2 H_\2 i m H 
m 2 ' Mvr^ 27 (16vr 2 ) 2 1 32 e 



-0.416(6) -i 1.000(7)]} , (3) 

where 7 = 0.577215664901532860607 is the Euler constant, and CI denotes 
the Clausen function, with C7(§) = 1.01494160640965362502. 

A number of checks were performed to make sure that these results are 
correct. 

The cancellation of logarithms and poles in £ was checked analytically for 
<5m^ 2 loop ^ a nd numerically for 8Zjj loop \ The complex integration paths 
for 5Z^ loop) were varied to make sure the integrand is indeed analytical 
along the integration path over the Feynman parameters (see ref. ||). The 
numerical value of SZ^ loop ^ agrees with the less precise one which can be 
derived from the momentum dependence of the Higgs selfenergy calculated 
in ref. ||. Note that the latter result implies numerical integrations over 
very different functions, and is sustained by the fact that its imaginary part 
agrees with the one-loop corrections to the Higgs width. The imaginary 
part of the momentum derivative of the Higgs selfenergy, which is given in 
eq. 3, 



- 1.000(7) (5 2 ^) 2 777^2 . < J -> 



obtained by numerically integrating the two-loop diagrams, agrees with the 
analytical result which can be derived by using the Cutkosky rule: 



2 m 2 H 2 1 3tt / Trx/3 5tt 2 \ 
{U m 2 w > (16vr 2 ) 2 4 ^ ^ 12 48 J 

= -1.002245142 (g 2 ^) 2 — ^ . (5) 
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Further checks were performed by calculating the two-loop selfenergy of 
the Goldstone bosons. After performing an expansion in the external mo- 
mentum, one can calculate the necessary Feynman diagrams analytically. 
The selfenergy of the Goldstone bosons at zero momentum agrees with the 
two-loop tadpole counterterm given in ref. ||, and the wave function renor- 
malization satisfies the generalized Ward identity at leading order in mjf, 
<5m^ = -m^SZ G , 

We are now in the position to calculate the two-loop radiative corrections 
to the Higgs fermionic width: 



Z H 


1 + bz%- loov) 


+ 5Z^- loop) 


+ ... 
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T" . ■ ■ , 



where 

Numerically, one obtains: 
The results are shown in fig. 2. 

Including the two-loop results reduces the magnitude of the radiative 
corrections, wich remain under 6% for the whole range of validity of the 
perturbative approach, up to ~ 1.4 TeV, where the two-loop correction 
exactely compensates the one-loop one. 

Finally, some comments on the possible relevance of perturbation the- 
ory beyond this limit are in order. In an attempt to sum the asymptotic 
series of eq. 6, fig. 3 shows the behaviour of the first nontrivial term of its 
Shanks transformation. Since the perturbative series is a power series, this 
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reduces actually to the Pade approximant [1/1]. Certainly, the convergence 
radius of the perturbative expansion is zero, and a diagonal sequence of Pade 
approximants might fail to sum it. To sum the perturbative series would 
require knowledge of its high order behaviour. On the much simpler case 
of the anharmonic oscillator one can show that finding the right nonlinear 
sequence transformation to sum an asymptotic series depends crucially on 
the knowledge of an estimate of the series' remainder fp~3|] . If such informa- 
tion is not available, the only thing one can do is to compare the results of 



different summation recipes [12 1 , but, of course, knowledge of only the one- 
loop and two-loop corrections does not leave much room for experimenting 
with nonlinear sequence transformations. 
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Figure captions 



Fig.l The topologies of the two-loop selfenergy diagrams. 

Fig. 2 The radiative corrections ^ to the fermionic Higgs width in the 
one-loop (solid line) and two-loop (dashed line) approximations as a func- 
tion of the mass of the Higgs boson. 

Fig. 3 The behaviour of the Pade approximant [1/1] (dotted line) in the 
strong coupling zone, compared to the one-loop (solid line) and two-loop 
(dashed line) approximations. 
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This figure "figl-l.png" is available in "png" format from: 
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